function [eta H Ac Bc] = RPfigures_a(pars)

global theta csi alpha delta s YK GK CK Xi
%Rea policy parameters
rho = pars(1);
gammaG = pars(2); gammaT = pars(3); lambdaG = pars(4); lambdaT = pars(5);

%Solving log-linear equations, can simplify resource constraint and EE 
%(substituting out for the interest rate and wages) to get:
% k' = V1*k + V2*c + V3*g + V4*tao
% E(c') = V5*c + V6*k' + V7*E(tao')
%with intermediate parameters:
V1 = 1 - delta + YK*(1-s)*alpha*(1+csi)/(alpha+csi);
V2 = - CK - YK*(1-s)*theta*(1-alpha)/(alpha+csi);
V3 = s*YK - GK;
V4 = -YK*(1-s)*(1-alpha)/(alpha+csi);
V5 = inv(1 + Xi*(1-alpha)/(alpha+csi));
V6 = - V5*Xi*(1-alpha)*csi/(alpha+csi);
V7 = -1/theta;
%Then using the policy rules for taxes and government spending get:
% k' = Z1*k + Z2*c + Z3*a
% E(c') = Z4*c + Z5*k' + Z6*a
Z1 = V1;
Z2 = V2 + V3*gammaG + V4*(gammaT + lambdaT*gammaG);
Z3 = (1+lambdaG)*(V3 + V4*lambdaT);
temp = 1 - V7*(gammaT + lambdaT*gammaG);
Z4 = V5*inv(temp);
Z5 = V6*inv(temp);
Z6 = V7*lambdaT*(1+lambdaG)*rho*inv(temp);
%But on this system, making guesses:
% k' = eta*k + H*a
% c = Ac*k + Bc*a
%matching coefficients, it is easy to reduce to a quadratic in eta:
% eta^2 - (Z1 + Z4 + Z2*Z5)*eta + Z1*Z4 = 0;
%which has roots:
eta1 = 0.5*(Z1 + Z4 + Z2*Z5 - sqrt(-4*Z1*Z4 + (Z1 + Z4 + Z2*Z5)^2));
eta2 = 0.5*(Z1 + Z4 + Z2*Z5 + sqrt(-4*Z1*Z4 + (Z1 + Z4 + Z2*Z5)^2));
%of which pick the smallest:
if abs(eta1)<abs(eta2); eta=eta1; else eta=eta2; end
%And plug this to solve tyhe other three equations for thwe other 3
% undetermined coefficients
Ac = (eta - Z1)/Z2;
Bc = (Z3*Z5 + Z6 - Ac*Z3)/(rho + Ac*Z2 - Z4 - Z2*Z5);
H = Z3 + Z2*Bc;
